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Single crystals of the dilute, rare earth bearing, pseudo-ternary series, Gda;Yi_2,Fe2Zn2o were 
grown out of Zn-rich solution. Measurements of magnetization, resistivity and heat capacity on 
Gda;Yi_a;Fe2Zn2o samples reveal ferromagnetic order of Gd^"*" local moments across virtually the 
whole series {x > 0.02). The magnetic properties of this series, including the ferromagnetic ordering, 
the reduced saturated moments at base temperature, the deviation of the susceptibilities from Curie- 
Weiss law and the anomalies in the resistivity, are understood within the frame work of dilute, S 
moments (Gd^^) embedded in a nearly ferromagnetic Fermi liquid (YFe2Zn2o). The s-d model is 
employed to further explain the variation of Tc with x as well as the temperature dependences of 
of the susceptibilities. 

PACS numbers: 75.10.Lp, 75.50.Cc, 75.20.Hr 
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I. INTRODUCTION 



Materials that are just under the Stoner limit mani- 
fest large electronic specific heat and enhanced paramag- 
netism and are sometimes known as nearly ferromagnetic 
Fermi liquids (NFFL)^'^. Archetypical examples, such 
as Pd^, NigGa^, TiBe^^ and YCoz^, have been studied 
for several decades. In addition to the interesting, in- 
trinsic properties of these compounds, the introduction 
of local moments into these highly polarizable hosts has 
lead to both experimental^ and theoretic interesli^i^. In 
such highly polarizable hosts, local moment impurities 
can manifest long range, ferromagnetic order even for 
very low concentrations (0.5 at.% Fe in Pd-" and 1 at.% 
Gd in Pd"). 

Recently, YFe2Zn2o was found to be a ternary exam- 
ple of a NFFL with a Stoner parameter Z ~ 0.9,— as 
compared to Z ~ 0.83 for Pd, indicating strongly cor- 
related electron behavior. When the large, S moment 
bearing, Gd^+ replaces the non-magnetic Y'^+ ions, it was 
found that GdFe2Zn2o has a remarkably high ferromag- 
netic Curie temperature(rc) of 86 K. Both of these com- 
pounds belong to the much larger, isostructural RT2Zn2o 
(R = rare earth, T = transition metal such as Fe, Co, 
Ni, Ru, Rh, Os, Ir, Pt)i^iii family, in which the R and 
T ions each occupy their own unique, single, crystallo- 
graphic sites. In these dilute, rare earth bearing, inter- 
metallic compounds (less than 5 at.% rare earth), the R 
ions are fully surrounded by Zn nearest and next nearest 
neighbors to form a Frank-Kasper-like Zn polyhedron; 
the T site is also surrounded by a nearest and next near- 
est neighbor, Zn shell. The shortest R-R spacing is ^6 
A. Motivated by these intriguing magnetic and structural 
properties, we focus, in this work, on the pseudo-ternary 
series Gda;Yi_j:Fe2Zn2o, which can be used as a model 
for studying the effects of titrating very dilute local mo- 
ments into a nearly ferromagnetic Fermi liquid. Given 



that RFe2Zn2o is a dilute, rare earth bearing intermetal- 
lic, dilution of Gd onto the Y site (i) changes the lattice 
parameter by less than 0.2 %, (ii) does not change the 
band filling, (iii) does not change the all Zn local envi- 
ronment of either the Gd or Fe ions, and (iv) allows for 
the dilution of Gd in the system to be studied down to 
X « 0.005, i.e. down to approximately 200 p.p. to. Gd. As 
shown below, single crystals of Gd2;Yi_a;Fe2Zn2o can be 
easily grown by a Zn, self flux metho d^^'^^ , and the Gd 
concentration can be consistently inferred via a variety 
of methods. 

In this paper, we report on the characterization of sin- 
gle crystals of Gda;Yi_2;Fe2Zn2o by X-ray diffraction. En- 
ergy Dispersive X-ray Spectroscopy (EDS), magnetiza- 
tion, resistivity and heat capacity measurements. These 
data reveal ferromagnetic order of the Gd^+ local mo- 
ment above 1.80 K for Gd concentration above x = 0.02. 
These results will be discussed within the framework of 
the so-called s-d model-^, based on the mean field approx- 
imation, and used to explain the variation of Tq across 
the series with respect to x. 



II. EXPERIMENTAL METHODS 

Single crystals of Gda;Yi_a;Fe2Zn2o were grown from a 
Zn-rich self flu x^^i^^ . For x > 0.02, high purity elements 
were combined in a molar ratio of (Gda;Yi_2:)2Fe4Zng4. 
For X less than 0.02, a Yp.gGdo.i master alloy was made 
via arc melting and appropriate amounts of this alloy 
were added to elemental Y to reduce the uncertainties as- 
sociated with weighing errors. The constituent elements 
(or alloy) were placed in a 2 ml, alumina crucible and 
sealed in a silica tube under approximately 1/3 atmo- 
sphere of high purity Ar (used to help reduce the evapo- 
ration and migration of zinc during the growth process) 
and then heated up to 1000 °C and cooled, over a period 
of 80 h, to 600 °C, at which point the remaining liquid 
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FIG. 1: Powder X-ray diffraction pattern of GdFe2Zn2o with 
a Si internal standard (using Cu Ka radiation) with main 
peaks indexed. Inset: the normahzed intensity of the (117) 
peak of Gdi:Yi_iFe2Zn2o for representative x values, with the 
positions calibrated by the nearby Si(002) peak. 



was decanted. Growths such as these often had only a 
few nucleation sites per crucible and yielded crystals with 
typical dimensions of at least 7x7x7 mm^ . Residual flux 
and/or oxide slag on the crystal surfaces was removed by 
using 0.5 vol.% HCl in H2O in an ultrasonic bath for 
1-2 h. The samples were characterized by room temper- 
ature powder X-ray diffraction measurements using Cu 
Ka radiation with Si (a = 5.43088 A) as an internal stan- 
dard in a Rigaku Miniflex powder diffractomcter(Fig. [1]). 
The Rietica, Rietveld refinement program was employed 
to obtain the lattice constants, which vary linearly for 
I < X < 0. This shift can be seen in the (117) peak 
position for selected x values (see Fig. [U inset). EDS 
measurements were made in a JEOL model 59101v-SEM 
with a Vantage EDS systerm for representative samples. 

In order to measure the electrical resistivity with a 
standard AC, four-probe technique, the samples were cut 
into bars using a wire saw. The bars typically had lengths 
of 2-4 mm parallel to the crystallographic [110] direc- 
tion, and widths and thicknesses between 0.2-0.4 mm. 
Electrical contact was made to these bars by using Epo- 
tek H20E silver epoxy, with typical contact resistances 
of about 1 Ohm. AC electrical resistivity measurements 
were performed with / = 16 Hz and I = 1 — 0.3 mA 
in a Quantum Design PPMS-14 or PPMS-9 instrument 
(T = 1.85 — 310 K). Temperature dependent specific heat 
measurements were also performed by using the heat 
capacity option of these Quantum Design instruments, 
sometimes using the '^He option. 

DC magnetization was measured in a Quantum Design 
superconducting quantum interference device (SQUID) 
magnetometer, in a variety of applied fields {H < 55 
kOe) and temperatures (1.85 K < T < 375 K). In some 
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FIG. 2: Magnetization M with respect to applied field H for 
a sample of Gdo.5Yo.5Fe2Zn2o at 150 K, 200 K and 300 K. 
The solid lines are guides to the eye. Inset a: detailed mag- 
netization of two samples of Gdo.5Yo.5Fe2Zn2o at 300 K. The 
data shown in black (same data as in main figure) has slope 
change feature (indicated by an arrow) ; while the data shown 
in red does not. Inset b: the difference of the black and red 
reveals the saturation of ferromagnetic impurity above 5 kOe. 



crystals, the magnetization with respect to magnetic field 
measurements at 300 K showed a slight non-linearity with 
a small slope change around 3 kOe (Fig. [2]). This spe- 
cific behavior is believed to be due to a small amount of 
ferromagnetic impurity, possibly Fe or FeOa; (2 x 10"^ 
/is/mol to 2 X 10^"^ /is/mol) on the crystal. This fea- 
ture is most likely extrinsic because the extent of the 
slope change is sample-dependent; some samples show- 
ing no feature at all (inset of Fig. [2]). This feature is 
most clearly seen when two samples from the same batch 
(one with feature, one without) are compared (Fig. [21 
inset a) or even subtracted from each other (Fig. [21 in- 
set b). Given that this small, extrinsic ferromagnetic 
contribution saturates by 7f « 10 kOe (Fig. [21 inset b), 
the high temperature susceptibility can be determined by 

ture region the intrinsic magnetization is a linear function 
of applied magnetic field for 20 kOe < H < 50 kOe (Fig. 
[2]). At lower temperatures, closer to Tq, the sample's in- 
trinsic magnetization become large enough that we can 
measure x{T) directly as M/H for H — 1 kOe. 



III. EXPERIMENTS RESULTS 

The size of the cubic unit cell, as determined by pow- 
der X-ray diffraction measurements, shows a linear de- 
pendence on X as it is varied from to 1 (Fig. [3]). The 
error bars of the lattice constants were estimated from 
the standard deviation determined by measurements on 
three samples from the same batch. These data are com- 
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FIG. 3: Gd concentration inferred from EDS (solid squares) 
and high temperature magnetic susceptibility (solid circles). 
The open triangles represent lattice constants. The dash line 
is location where inferred x equals nominal x and also repre- 
sents a linear dependance of the lattice parameter. 



pliant with Vegard's law and imply that the nominal x 
is probably close to the actual x. 

In order to check this further, EDS was used. This is a 
direct method to determine x, although it loses some of 
its accuracy because of the low, total rare earth concen- 
tration (< 5 at.%). Nevertheless, several representative 
members of the Gda;Yi„a;Fe2Zn2o series were measured 
and the inferred x values are close to the nominal x values 
within the fairly large error bars (Fig. [3]). 

Another way to estimate the concentration of gadolin- 
ium in the grown crystals is based on the analysis of the 
high temperature magnetic susceptibility data, which can 
be expressed as: 



(1) 



Experimentally, XGd'-^+ obeys the Curie- Weiss law above 
150 K (Fig. 2^), from which the paramagnetic Curie 
temperature 6c and Curie constants C can be extracted. 
The value of x can be inferred by fixing the effective 
moment of Gd^+ as 7.94 /i^. These values of inferred x 
are also plotted in Fig. [S] The agreement between each 
of these three different methods of determining inferred 
X and the nominal x value is good and for the rest of this 
paper nominal values will be used to estimate actual Gd 
content. 

Another aspect of Fig. 0] that is noteworthy is that 
all XGd^+ data sets deviate from their high tempera- 
ture Curie- Weiss behaviors as the system approaches the 
magnetic ordering temperature. Since high fields can 
shift and broaden the features associated with ferromag- 
netism, at lower temperatures a field of 1 kOe was used 
(Fig. [3|d). Whereas this deviation cannot be associated 
with the formation of superparamagnetic clusters above 
Tq (this would cause a slope change toward the hori- 
zontal rather than toward the vertical), it can be un- 
derstood in terms of an increasing coupling between the 
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FIG. 4: l/xGd3+ vesus temperature for representative mem- 
bers of the Gda;Yi_a;Fe2Zn2o sereies. Note: data is normalized 
to mole Gd using x inferred from high-temperature data, (a) : 
obtained under high magnetic field, (b) Solid lines: obtained 
under 1 kOe applied field; dash lines: under high magnetic 
field. 



Gd^^^ local moments associated with the strongly tem- 
perature dependent, polarizable electronic background of 
the YFe2Zn2o matri3slS.(see discussion below). 

Figure [5] shows the temperature dependent magneti- 
zation in an external field H = 1000 Oe for the whole 
range of x values. Ferromagnetic ordering can be clearly 
seen below 90 K for x = 1. This ordering tempera- 
ture decreases monotonically as x decreases, although 
the exact values of Tq can not be unambiguously inferred 
from these plots. For x < 0.035, it becomes difficult to 
determine whether the compounds manifest ferromag- 
netism above the base temperature (1.85 K) based on 
the M{T) curves alone. Even at 1000 Oe, for x > 0.25, 
the low-temperature magnetization is just slightly below 
the Hunds ground state value 7 /is/Gd at the base tem- 
perature (Fig. [5^). For x < 0.25 the low temperature, 
H — 1000 Oe, magnetization decreases with decreasing 
X (Fig. [Hd). 

Field-dependent magnetization measurements were 
made for each sample at base temperature (Fig. [6]). For 
compounds with x > 0.035, the magnetization rapidly 
saturates as the magnetic field increases, consistent with 
a ferromagnetic ground state at 1.85 K. For x < 0.01, 
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FIG. 5: Temperature dependent magnetization of 
Gd^Yi-^Fe2Zn2o, H = 1000 Oe, for (a) 1.0 > x > 0.175, (b) 
X < 0.175. 



the M{H) curves vary more smoothly with H and are 
more consistent with a paramagnetic state at 1.85 K. 
The X — 0.02 data are more ambiguous and require a 
still more detailed analysis. 

For H > IQ kOe the M{H) data for x < 0.05 vary 
approximately linearly with H and have slopes compa- 
rable to that of YFe2Zn2o.(Fig. [6)3) For all x values the 
magnetization can be thought of as a combination of the 
magnetization of Gd'^"'" ions and the highly polarizable 
background. In order to extract the magnetization of the 
Gd^+ ions, a background of MYFe2Zn2o was subtracted 
from the M{H) data. The MaiH) data are plotted in 
Fig. [6]: normalized to the nominal x values. For x > 0.25 
the saturated magnetization is essentially constant with 
a value slightly less than 7 /is/Gd.^^ For x < 0.25 there 
is an apparent decrease in the saturated magnetization 
with decreasing x, but it should noted that the error 
bars, coming from the estimated ±0.02 uncertainty of x, 
increase with decreasing x. These increasing error bars 
make it unclear whether the saturated moment of the Gd 
impurities is constant or decreasing in the small x limit. 

A fuller analysis of M{H) data, particularly the anal- 
ysis of magnetization isotherms known as Arrott plots^^, 
at a set of temperatures near Tq has been found to be a 
useful, and in some case even the best method to deter- 
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FIG. 6: (a) and (b) Field dependent magnetization of 
Gd2,Yi_a;Fe2Zn2o at 1.85 K. (c) Field dependent magnetiza- 
tion of Gda;Yi_a;Fe2Zn2o at 1.85 K, normalize to Gd"^""" con- 
tent (see text). The error bars were estimated by allowing for 
a ±0.02 variation of x 



mine Tc for the x < 0.25 samples. The method is based 
on the mean field theory, in which is linear in //M 
with zero intercept at the critical temperature Tc , where 
/ is the internal field, equal to the difference between the 
external, applied field H and the demagnetizing field Dm- 
For an ellipsoid of Gd:EYi_2;Fe2Zn2o, the demagnetizing 
field equalsi^: 
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47r 



Na = 0.061 ■ D- M (2) 



where M is the magnetization (emu/mol); Z? is a geomet- 
ric factor that can range from 1 to 0, and is Avogadro 
number. Thus I/M, in units of /i^/kOe , is: 
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Using H, instead of / , in Arrott plots will shift the data 
along H/M axis in the positive direction by 0.34 x D, 
which would experimentally introduce an error in the 
value of Tc for a flat shaped sample {D ^ 1) of 
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FIG. 7: Arrott plots for representative members of the 
Gd:„Yi_^Fe2Zn2o series: x = (a) 0.5, (b) 0.035 and (c) 0.02. 



FIG. 8: Zero-field resistivity for current along the [110] direc- 
tion. The arrows represent Tc determined from Arrott plot 
analyses. 



GdFe2Zn2o. Nevertheless, even in this extreme case, this 
error drops as x decreases due to reduction of the sam- 
ples' magnetization as Gd'^+ is diluted out (notice the 
different scale of the M axis for x < 0.05 in Fig. [71). Due 
to these concerns, rod-like-shape samples were measured 
along their long axis for the magnetization isotherms for 
samples with x > 0.5. This shape ensures D is min- 
imized. Figures [7^ and b show Tc = 57 ± 0.5 K for 
X = 0.5 and Tc = 4.5 ± 0.5 K for x = 0.035 respec- 
tively. For X = 0.02, Fig. [7t shows Tc = 1.85 K, a result 
that helps explain the difficulty experienced in determin- 
ing the base-temperature magnetic state based on the 
M{T) and M{H) data discussed above. The Tc values 
determined for the Arrot plot analysis for all x are shown 
below in Fig. [T3l 

The temperature dependent electric resistivity data, 
p{T) (measured in zero applied magnetic field), of the 
Gd2;Yi_a;Fe2Zn2o compounds are shown, for representa- 
tive X values, in Fig. [51 For x > 0.25, p{T) curves show 
a kink at Tc due to the loss of spin disorder scattering 
below this temperature. In contrast, for x < 0.175, no 
clear kink can be detected. Tc values deduced from the 
maximum of dp/dT (not shown here) are compatible with 
the values obtained from the Arrott plots (see Fig. [T5b 
below) . 

Further information can be extracted from the the 
Gd3;Yi_2;Fe2Zn2o p(T) data by assuming that the total 



resistivity of the compound can be written as: 

P{T) = Po + Pph{T) + Pmag{T), (4) 

where po is a temperature independent, impurity scatter- 
ing term, pph is the scattering from phonons and Pmag is 
scattering associated with the interaction between con- 
duction electrons and magnetic degrees of freedom. In 
this series of pseudo-ternary compounds, the high tem- 
perature (T ^ Tc) phonon contribution, pph, should be 
essentially invariant (due to the very dilute nature of the 
R ions). The magnetic contribution to the resistivity, 
Pmag i wiU be the combination of contributions from con- 
duction electron scattered by (i) the 4f local moments 
and (ii) the spin fluctuations of 3d electrons (from Fe 
sites), both of which should saturate in the high tem- 
perature limit. Based on the analysis above, the high 
temperature resistivity of the whole series should be sim- 
ilar (modulo an offset) and manifest similar slopes due to 
the electron-phonon scattering. This is indeed the case: 
the data show linearity of p(T) above 250 K with the 
slopes differing by less than 8%; less than the estimated 
dimension error (10%) of these bar-like-shape samples. 

The magnetic and disorder contribution to the resistiv- 
ity can be estimated by (i) removing the geometric error 
by normalizing the high temperature slope of all p(T) 
plots to that of YFe2Zn2o and then (ii) subtracting the 
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Py{T) data from the p normalized data. 
The normahzed p is given as: 

and 

^P ~ PGd^normalized — PY ■ (6) 

The resulting Ap will not only show the conduction 
electron scattering from the 4f local moments, but will 
also include scattering associated with the interaction 
between the 4f local moment and 3d electrons, espe- 
cially near Tq. The temperature dependent Ap curves 
for the Gda;Yi_a;Fe2Zn2o compounds are presented in 
Fig. [21 A pronounced upward cusp is centered about 
Tc for X > 0.25. For x < 0.25 the loss of spin disor- 
der feature becomes harder (or even impossible) to re- 
solve, but the enhanced scattering above Tc persists. 
The decrease of Ap with T below Tq is a common in 
ferromagnetic systems and can be explained as the result 
of a loss of spin disorder scattering of conduction elec- 
trons. On the other hand, the behavior of Ap above Tc 
must come from a different conduction electron scatter- 
ing process. A similar feature in Ap is found in RFe2Zn2o 
(R = Tb-Er) for T > To^, but not in isostructural 
GdCo2Zn2o which orders antiferromagnetically at a much 
lower temperature^^. 

The specific heat of the Gd2.Yi_a;Fe2Zn2o compounds 
(Fig. [TO]) can be thought of as the sum of the contribu- 
tions from electronic, vibrational and magnetic degrees of 
freedom. To remove the vibrational and electronic parts 
(at least approximately), the specific heat of YFe2Zn2o 
and LuFe2Zn2o were used to estimate the background. 
The assumption that YFe2Zn2o and LuFe2Zn2o closely 
approx the non-magnetic Cp of the Gd2;Yi_a;Fe2Zn2o 
series is supported by the fact that the difference be- 
tween the measured Cp of YFe2Zn2o, LuFe2Zn2o and 
Gd2,Yi_a;Fe2Zn2o in the temperature region 20 K higher 
than Tc is on the order of a percent. Since LuFe2Zn2o has 
a molar mass closer to that of GdFe2 Zn2o than YFe2 Zn2o , 

the combination of {x)CLuFe2Zn2o + (1 - x)CYFe2Zn2o 

is thought to be even closer to the non-magnetic back- 
ground of CGd^Yi^^Fe2Zn2o- 

Figure [11] shows 

AC = CGd^Yi-^Fe2Zn2o i"^) 
— {x)C LuF e2Zn2o " (1 ~ x)CYFe2Zn2o 

for X > 0.175 (a) and < 0.175 (b), where the arrows indi- 
cate the Tc values determined from the Arrott plot anal- 
yses. The magnetic ordering manifests itself as a broad 
feature in AC with Tc occurring at, or near, the max- 
imum slope. Figure [TO] shows that this feature persists, 
relatively unchanged in shape, down to a; = 0.1. For val- 
ues of X < 0.1 the feature broadens further, but is still 
distinct. This shape of AC is not unusual for Gd-base 
intermetallics with ferromagnetic order; for example, a 
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FIG. 9: Temperature variation of Ap (see text). The arrows 
represent Tc determined from Arrott plot analysis of magne- 
tization measurements. 



similar feature is seen in GdPtIn (Tc ~ 68 K )^. It 
should be noted that this AC feature is distinct from 
that associated with a spin-glass freezing: the maxima 
all occur at or below Tc, whereas a spin glass manifests 
a broad peak above the freezing temperature^^ . 

The X dependence of the paramagnetic Curie tempera- 
ture (Oc), ferromagnetic ordering temperature (Tc) and 
saturated moments per Gd (p-sat) for each x are shown 
in Fig. [TOk . b and c respectively. The values of the mag- 
netic entropy, estimated by Sm = / ^^dT , are shown in 
Fig. [TOH . Both Oc and Tc decrease monotonically with 
X. At first glance, the negative values of 9c for x < 0.25 
are unexpected and seem to be in contradiction with the 
existence of ferromagnetic ground state. However, the 
data analyses in Eq. [1] ignores, the increasingly strong, 
polarizable background associated with the near Stoner 
limit conduction electrons at intermediate temperatures. 
Furthermore, as shown in Fig. [Dd, this low temperature 
effect becomes even more pronounced for small x. Al- 
though, as discussed earlier, the uncertainty of x makes 
the x-variation of psat ambiguous for small x, even the 
large x members of the Gd2:Yi_a;Fe2Zn2o series manifest 
reduced saturated moments. This is attributed to the in- 
duced moment on the 3d electrons, which is anti-parallel 
to the Gd moment The magnetic entropy shown in 
Fig. [TOh associated with the ordered state is equal to. 
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FIG. 10: Temperature variation of specific lieat Cp of the 
Gd2:Yi_a;Fe2Zn2o series for x = 0.75, 0.5 and 0. The arrows 
represent Tc determined from Arrott plot analyses. 



or slightly larger than, the magnetic entropy associated 
with the Hunds ground state of Gd3+(S = 7/2). This 
fact indicates that the main part of the magnetic specific 
heat of the series of \-x^%2ii^iQ is the contribution 

from the magnetic degrees of freedom of the Gd'^"'" lo- 
cal moments. The contribution to the magnetic specific 
heat from the itinerant electrons probably exists, but is, 
at most, comparable with the measurement uncertainty. 
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IV. ANALYSIS AND DISCUSSION 

For rare earth bearing intermetallics, the interaction 
between 4f local moments is primarily mediated by means 
of polarization of the conduction electrons. Regard- 
less of the details of the mechanism involved in this 
interactio n^^'^'^ , we emphasis that the 3d electrons from 
Fe sites act as important mediators of the Gd-Gd interac- 
tion in Gd2;Yi_a;Fe2Zn2o system. In YFe2Zn2o, the inter- 
action between 3d electrons is not sufficient to split the 
conduction band but is large enough to make the com- 
pound exhibit strongly enhanced paramagnetism. When 
Y'^^ ions are fully replaced by Gd"^"*" ions, these 3d elec- 
trons are polarized by the Gd'^"'" local moments. The 
interaction between 3d electrons assists in stabilizing the 
splitting of the conduction electron band and enhances 
the magnetic interaction between Gd"^"*" local moments, 
resulting in the remarkably high, ferromagnetic transi- 
tion temperature for GdFe2Zn2o. This physical picture 
is consistent with the results of the band structure calcu- 
lation which predicts the Fe induced moment as 0.67/^^ 
in the ground state of GdFe2Zn2cr^ii^. 

In order to perform further analysis on the magnetic 
properties of Gd2:Yi_a;Fe2Zn2o, a comparision with the 
binary RC02 (R = rare earth) intermetallics is useful. 
YC02 and LuCo2 show nearly ferromagnetic behavior 



while the series of compounds, (Gd-Tm)Co2, with 4f local 
moments manifest a ferromagnetic ground stat o^'*'^^ . In 
addition to these magnetic similarities, the resemblance 
between the crystal structure of RT2Zn2o and the so- 
called C-15 Laves structure of RCo2^is noticeable: both 
rare earth and transition metal ions occupy same unique, 
single crystallographic sites in same space group: FdSm. 
The unit cell of the RT2Zn2o compounds can be thought 
of as an expansion of the C-15 Laves phase unit cell via 
the addition of a large number of Zn ions. 

Well-studied for several decades, the series of (Gd- 
Tm)Co2 has been treated as an example of 4f local mo- 
ments embedded in a nearly ferromagnetic host: YC02 
or LuCo2. The so-called s-d model has been employed 
by Bloch and Lemaire^^ and Bloch et. ali^ to explain 
their magnetic properties. This model was first intro- 
duced by Takahashi and Shimizu^^ to understand the 
magnetic properties of alloys of the nearly ferromagnetic 
transition metal, Pd, with dilute Fe or Co local moment 
impurities. In this model, the polarization effect of the 
local moments on the itinerant electrons is considered in 
terms of a molecular field. Motivated by the similarity 
of the magnetic properties and the crystal structure of 
RFe2Zn2o and RC02, we applied the s-d model to the 
Gd2;Yi_2.Fe2Zn2o series. 
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This model considers one magnetic system consist- 
ing of two types of spins: one local moment, and the 
other giving rise to an exchange-enhanced paramagnetic 
susceptibility.^ - Assuming the interaction between Gd lo- 
cal moments is only via the conduction electrons, we ap- 
ply this model to Gd2;Yi_a;Fe2Zn2o- Under an applied 
field H, for T > Tc, the magnetization of the Gd local 
moments and the conduction electrons are: 



Mad = {xCGd/T){H + UGd-eM,) 



Xefi{H + Ue-eMe + ncd-eMod) 



(8) 
(9) 



where Ccd is the Curie constant of the Gd'^+ local mo- 
ments; ncd-e, n^-e are molecular-field coefficient repre- 
senting the interaction between itinerant electrons and 
Gd^+ local moments, and itinerant electrons with them- 
selves, respectively; Xe,o is the paramagnetic susceptibil- 
ity without exchange enhancement. The total magneti- 
zation of Gd2;Yi_a;Fe2Zn2o is the sum of Mcd and Me. It 
should be noted that when a; = 0, the total susceptibility 
reduces to the exchange-enhanced susceptibility: 



Xe = XYFe-iZn-io 



Me 

H 



1 — n, 



e — eXe,0 



(10) 



which is simply the Stoner enhanced susceptibility of 
YFe2Zn2o. 



FIG. 13: (a) Paramagnetic Curie temperature, 9c, (b) Fer- 
romagnetic ordering temperature, Tc, (c) saturated moment 
per Gd, j-isat and (d) magnetic entropy, Sm with respect to 
X for the Gd2,Yi_a;Fe2Zn2o series. The values of Tc in (b) 
were determined by Arrott plot analyses (black circle) and 
the resistivity measurements(red circle). The solid line in (d) 
represents Sm = xRlnS (R is gas constant), the magnetic 
entropy of Gd^^ Hunds ground state. The error bars are es- 
timated as 1 % of the total entropy, 5* = J^'~' 



Assuming that the electronic structure of the con- 
duction band and the position of the Fermi level in 
the paramagnetic state are the same across the whole 
Gda;Yi_a;Fe2Zn2o series, from Eqs. [5HTU1 one gets the 
total susceptibility of Gda;Yi_a;Fe2Zn2o 



xCcd 



T - XY Fe2Zn2o'"-Gd~e^^Gd 
XYFe2Zn2o{T + 2nGd-eXCGd) 



(11) 



T - XYFe2Zn2a'nGd- 



.xCGd 



If one assumes the coupling between the pure spin mo- 
ment (S ~ 7/2) of the Gd'^^ and the conduction electron 
spin cr ((T = 1/2) to be a Heisenberg exchange interac- 
tion, 2 Jos' ■ CT, where Jo is the exchange parameter, then 
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FIG. 14: Tc of Gd:j,Yi_^Fe2Zn2o versus x ■ XYFe2Zn2o{Tc)- 
The solid line is linear fit through the origin point 
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FIG. 15: 1/x of GdiYi_a;Fe2Zn2o versus T for representative 
X values. Dot lines: measured under 1 kOe applied filed; 
dash lines: obtained under high magnetic field; solid lines: 
calculated results. (See text) 



the molecular field coefficient 



ncd-e 



^Jo/(2m|7V) 



(12) 



where N is the number of the rare earth ions per volume. 

The Gd2;Yi_2;Fe2Zn2o system will become ferromag- 
netic when XGd^Yi^^Fe2Zn2o diverges. Thus, 



Tc — XYFe2Zn2o{Tc)n%cL~e^CGd 



(13) 



^ ' XY Fe2Zn2o (Tc) 



J§SiS+l) 



where fc^ is the Boltzmann constant. 

Equation [13] reveals that Tq depends on the product 
of x and XYFe2Zn2o{Tc)i rather than just x. This is con- 
sistent with Fig. 113b showing a nonlinear dependence of 
Tc on X. Figure [Til shows that the values of Tc depend 
linearly on the product xxy Fe2Zn2o{Tc) across the whole 
series. From Fig. [Hthe slope equals 2.955 ±0.0037 x lO"* 
K mol/emu and thus Jq can be extracted as 3.96 ± 0.05 
meV. 

In addition to the magnetic ordering, this model can 
also explain the curious temperature dependence of the 
1/x versus T data for the Gd2,Yi_a;Fe2Zn2o series. Set- 
ting Jo = 3.96 meV, one obtains the temperature de- 
pendent, total susceptibility of Gda;Yi_2,Fe2Zn2o. The 
results of l/xGd^Yi^^Fe2Zn2a for representative x values 
are shown as the solid lines in Fig. [T51 whereas the dotted 
lines and the dash lines present the experimental results 
under 1 kOe and high magnetic field, representatively. 
These calculated results qualitatively reproduce the ex- 
perimental, temperature dependent susceptibilities, es- 
pecially their deviation from the Curie- Weiss law close 
to Tc. It should be noted that the x data in Fig. [15] 
is the full X without any subtraction of " non-magnetic" 
background. In this sense Fig. [TS] and the s-d model, 



appear to treat the magnetization data more fully than 
the simple assumption behind Eqn. [1] 

In addition to the thermodynamic properties discussed 
above, the feature in Ap above Tc (Fig. [9]) is also worth 
discussing further. The upward-pointing cusp at Tc of 
Ap(T) is associated with the sign change of d/S.p/dT , 
from negative to positive as the temperature decreases. 
This feature is absent from simple models of p(r) 30,3i ^ 
based on the models assuming a single lattice of mag- 
netic ions and a single band of conduction electrons. This 
theoretical model is over-simplified for Gd3;Yi_a;Fe2Zn2o, 
a strongly correlated electron system. Similar unusual 
upward cusps in Ap(T) at Tc were found in the elec- 
tric transport measurements of RCo2'^^. They were ex- 
plained by invoking an increasing, non-uniform fiuctuat- 
ing f-d exchange interaction, which provides an increase 
of spin fluctuation of 3d-electron subsystem as the tem- 
perature approaches Tc in paramagnetic state, which 
in turn leads to increased conduction electron scatter- 
ing. The Gda;Yi_a;Fe2Zn2o system (and indeed the other 
RFe2Zn2o compoundsi^) present another, clear example 
of this behavior. 



V. SUMMARY 

We presented a set of data including magnetization, 
electrical transport and specific heat, measured on flux- 
grown single crystals of Gda;Yi_a;Fe2Zn2o. We found that 
the series order ferromagnetically above 1.85 K for x > 
0.02. The variation of Tc with respect to x, as well as the 
curious temperature dependent magnetic susceptibilities, 
are well explained by a modifucation of the s-d model 
based on molecular field approximation. The enhanced 



10 



electronic scattering found for T > Tc is qualitatively 
consistent with this model and similar to that found for 
the related RC02 compounds. 
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